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Abstract 

We discuss the current status of theoretical and experimental constraints on the real Higgs singlet extension 
of the Standard Model. For the second neutral (non-standard) Higgs boson we consider the full mass 
range from 1 GeV to 1 TeV accessible at past and current collider experiments. We separately discuss 
three scenarios, namely, the case where the second Higgs boson is lighter than, approximately equal to, or 
heavier than the discovered Higgs state at around 125 GeV. We investigate the impact of constraints from 
perturbative unitarity, electroweak precision data with a special focus on higher order contributions to the 
W boson mass, perturbativity of the couplings as well as vacuum stability. The latter two are tested up 
to a scale of ~ 4 x 10 10 GeV using renormalization group equations. Direct collider constraints from Higgs 
signal rate measurements at the LHC and 95% C.L. exclusion limits from Higgs searches at LEP, Tevatron 
and LHC are included via the public codes HiggsSignals and HiggsBounds, respectively. We identify the 
strongest constraints in the different regions of parameter space. We comment on the collider phenomenology 
of the remaining viable parameter space and the prospects for a future discovery or exclusion at the LHC. 
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I Introduction 


I. INTRODUCTION 

The LHC discovery mm of a Higgs boson in July 2012 has been a major breakthrough in modern 
particle physics. The first runs of the LHC at 7 and 8 TeV are now completed and the main results 
from various experimental analyses of the Higgs boson properties have been presented at the 2014 
summer conferences. So far, the discovered state is well compatible isffia with the interpretation 
in terms of the scalar boson of the Standard Model (SM) Higgs mechanism |11H15] . A simple 
combination of the Higgs mass measurements performed by ATLAS [16] and CMS m yields a 
central value of 


m H = (125.14 ± 0.24) GeV. (1) 

If the discovered particle is indeed the Higgs boson predicted by the SM, its mass constitutes 
the last unknown ingredient to this model, as all other properties of the electroweak sector then 
follow directly from theory. The current and future challenge for the theoretical and experimental 
community is to thoroughly investigate the Higgs boson’s properties in order to identify whether 
the SM Higgs sector is indeed complete, or instead, the structure of a more involved Higgs sector is 
realized. On the experimental side, this requires detailed and accurate measurements of its coupling 
strengths and CV structure at the LHC and ultimately at future experimental facilities for Higgs 
boson precision studies, such as the International Linear Collider (ILC) [18] . A complementary 
and equally important strategy is to perform collider searches for additional Higgs bosons. Such 
a finding would provide clear evidence for a non-minimal Higgs sector. This road needs to be 
continued within the full mass range that is accessible to current and future experiments. 

In this work, we consider the simplest extension of the SM Higgs sector, where an additional real 
singlet field is added, which is neutral under all quantum numbers of the SM gauge group [19, 20] 
and acquires a vacuum expectation value (VEV). This model has been widely studied in the 
literature f2TH48j . Here, we present a complete exploration of the model parameter space in the 
light of the latest experimental and theoretical constraints. We consider masses of the second 
(non-standard) Higgs boson in the whole mass range up to 1 TeV, thus extending and updating 
the findings of previous work [41]. This minimal setup can be interpreted as a limiting case for 
more generic BSM scenarios, e.g. models with additional gauge sectors [49] or additional matter 
content m Ell¬ 
in our analysis, we study the implications of various constraints: We take into account bounds 
from perturbative unitarity and electroweak (EW) precision measurements, in particular focussing 
on higher order corrections to the W boson mass [43] . Furthermore, we study the impact of 
requiring perturbativity, vacuum stability and correct minimization of the model up to a high 
energy scale using renormalization group evolved couplings 1 . We include the exclusion limits from 
Higgs searches at the LEP, Tevatron and LHC experiments via the public tool HiggsBounds [52]- 


The value of this high energy scale is chosen to be larger than the energy scale where the running SM Higgs quartic 
coupling turns negative. This will be made more precise in Section 


III 
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II The model 


A Potential and couplings 


[55] . and use the program HiggsSignals [56] (cf. also Ref. [57]) to test the compatibility of the 
model with the signal strength measurements of the discovered Higgs state. 

We separate the discussion of the parameter space into three different mass regions: (i) the 
high mass region, mu £ [130,1000] GeV, where the lighter Higgs boson h is interpreted as the 
discovered Higgs state; (ii) the intermediate mass region, where both Higgs bosons h and H are 
located in the mass region [120,130] GeV and potentially contribute to the measured signal rates 
and (in) the low mass region, G [1,120] GeV, where the heavier Higgs boson H is interpreted 
as the discovered Higgs state. 

We find that the most severe constraints in the whole parameter space for the second Higgs mass 
m H < 300 GeV are mostly given by limits from collider searches for a SM Higgs boson as well as by 
the LHC Higgs boson signal strength measurements. For mu > 300 GeV limits from higher order 


contributions to the W boson mass prevail, followed by the requirement of perturbativity of the 
couplings which is tested via renormalization group equation (RGE) evolution. For the remaining 
viable parameter space we present predictions for signal cross sections of the yet undiscovered 
second Higgs boson for the LHC at center-of-mass (CM) energies of 8 and 14TeV, discussing both 
the SM Higgs decay signatures and the novel Higgs-to-Higgs decay mode H —»• hh. We furthermore 
present our results in terms of a global suppression factor k for SM-like channels as well as the 
total width T of the second Higgs boson, and show regions which are allowed in the (/c,r) plane. 

The paper is organized as follows: In Section [TT] we briefly review the model and the chosen 
parametrization. In Section III we elaborate upon the various theoretical and experimental con¬ 
straints and discuss their impact on the model parameter space. In Section IV a scan of the full 
model parameter space is presented, in which all relevant constraints are combined. This is followed 
by a discussion of the collider phenomenology of the viable parameter space. We summarize and 
conclude in Section IVl 


II. THE MODEL 

A. Potential and couplings 

The real Higgs singlet extension of the SM is described in detail in Refs. praSDES]. Here, we 
only briefly review the theoretical setup as well as the main features relevant to the work presented 
here. 

We consider the extension of the SM electroweak sector containing a complex SU(2)l doublet, 
in the following denoted by <h, by an additional real scalar S which is a singlet under all SM gauge 
groups. The most generic renormalizable Lagrangian is then given by 

D+ d^Sdf.S - V($, S ), (2) 
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II The model 


A Potential and couplings 


with the scalar potential 


V($,S) 



= -m 2 ^ - n 2 S 2 + AiC^^) 2 + A 2 S 4 + X 3 &$S 2 . 


(3) 


Here, we implicitly impose a Z 2 symmetry which forbids all linear or cubic terms of the singlet held 
S in the potential. The scalar potential V (<f>, S) is bounded from below if the following conditions 
are fulfilled: 


4AiA 2 — A 2 > 0, 
Ai,A 2 > 0, 


(4) 

(5) 


cf. Appendix |A} If the first condition, Eq. Q, is fulfilled, the extremum is a local minimum. The 
second condition, Eq. ([5]), guarantees that the potential is bounded from below for large held 
values. We assume that both Higgs helds <f> and S have a non-zero vacuum expectation value 
(VEV), denoted by v and x, respectively. In the unitary gauge, the Higgs helds are given by 


$ = 



5 = 


V V2 


h! + x 


Expansion around the minimum leads to the squared mass matrix 

2Aiu 2 X 3 vx \ 


M 2 = 


X 3 vx 2A 2 .t 2 


with the mass eigenvalues 


ml = Aiu 2 + A 2 x 2 - \J (Aiu 2 - A 2 x 2 ) 2 + (A 3 xu) 2 , 
rn 2 H = Xiv 2 + A 2 x 2 + y 7 (Aiu 2 - A 2 x 2 ) 2 + (X 3 xv) 2 , 


( 6 ) 

(7) 

( 8 ) 

(9) 


where h and H are the scalar helds with masses m/j and mu respectively, and m? < m? H by 
convention. Note that > 0 follows from Eq. @ and we assume Eqs. 0 and ([5]) to be fulfilled 
in all following definitions. The gauge and mass eigenstates are related via the mixing matrix 



cos a — sin a 
sin a cos a 



where the mixing angle — § < a < f is given by 


sin 2 a = 

cos 2a = 


X 3 xv 

V ' (Aiu 2 - A 2 x 2 ) 2 + {X 3 xv) 2 ’ 

A 2 .t 2 — Aiu 2 

1 /(Aiu 2 - X 2 X 2 ) 2 + (X 3 xv) 2 


( 10 ) 


( 11 ) 

( 12 ) 
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A Potential and couplings 


It follows from Eq. (10) that the light (heavy) Higgs boson couplings to SM particles are suppressed 
by cos a (sin a). 

Using Eqs. <§ and we can express the couplings in terms of the mixing angle a, the 
Higgs YEVs x and v and the two Higgs boson masses, m,h and mu'- 


1 m h , (' m H - m l) ...2 
Al = 2i? + 


mr 


m 


H 


2v 2 

,™ 2 (m 2 H -m 2 ) 

X2 ~ 2^ + 2xA 


sin" a = —ty cos a H-sin" a, 

2v z 2v z 


mi 


m 


H 


cos“ a = — ,J k sm a H-^ cos“ a, 

2x z 2x z 


A 3 = 


(m 


H 


~ m l ) 


2vx 


sin (2a). 


(13) 


If kinematically allowed, the additional decay channel H —>• hh is present. Its partial decay 
width is given by 


Iff- 


>hh — 




/12 


8irm,H 


T- 


4m? 


m 


2 > 

H 


where the coupling strength p' of the H —> hh decay reads 


h = 


sin (2a) 
2vx 


m 


(sin av + cos a x) ( rn^ H—) . 


(14) 


(15) 


We therefore obtain as branching ratios for the heavy Higgs mass eigenstate 

^ H^hh 


HR h- 


>hh — 


tot 


. 2 TsM./f-s-SM 
HR//—ySM = sm a x -—-, 

-L tot 


(16) 


where TsM.ff^SM denotes the partial decay width of the SM Higgs boson and H —> SM represents 
any SM Higgs decay mode. The total width is given by 

Ttot = sin 2 a x Tsm, tot + ^H^hh, 

where Tsm, tot denotes the total width of the SM Higgs boson with mass mu- The suppression by 
sin 2 a directly follows from the suppression of all SM-like couplings, cf. Eq. (10). For // = 0, we 
recover the SM Higgs boson branching ratios. 

For collider phenomenology, two features are important: 


• the suppression of the production cross section of the two Higgs states induced by the mixing, 
which is given by sin 2 a (cos 2 a) for the heavy (light) Higgs, respectively; 

• the suppression of the Higgs decay modes to SM particles, which is realized if the competing 
decay mode H —>• hh is kinematically accessible. 


For the high mass (low mass) scenario, i.e. the case where the light (heavy) Higgs boson is iden¬ 
tified with the discovered Higgs state at ~ 125 GeV, |sina| = 0(1) corresponds to the complete 
decoupling of the second Higgs boson and therefore the SM-like scenario. 
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B. Model parameters 


At the Lagrangian level, the model has five free parameters, 

Ai, A 2 , A 3 , v, x, 


(17) 


while the values of the additional parameters /i 2 , m 2 are fixed by the minimization conditions to 
the values 


m = Ai v 2 + 


2 1 A3 2 
-ar, 


2 \ 2 1 A3 2 

/r = A 2 x + —v , 


(18) 

(19) 


cf. Appendix [Aj In this work, we choose to parametrize the model in terms of the independent 
physical quantities 


777 ,/^, mu, a, v , tan/3 = —. 

x 


( 20 ) 


The couplings Ai,A 2 and A 3 can then be expressed via Eq. (13). The vacuum expectation value 
of the Higgs doublet is given by the SM value v ~ 246 GeV. Unless otherwise stated, we 
fix one of the Higgs masses to be = 125.14 GeV, hence interpreting the Higgs boson h/H 

as the discovered Higgs state at the LHC. In this case, we are left with only three independent 
parameters, m = mjj/h, sin a, tan/3, where the latter enters the collider phenomenology only via 
the additional decay channel 2 H —>• hh. 


III. THEORETICAL AND EXPERIMENTAL CONSTRAINTS 


We now discuss the various theoretical and experimental constraints on the singlet extension model. 
In our analysis, we impose the following constraints: 

( 1 .) limits from perturbative unitarity, 

(2.) limits from from EW precision data in form of the S, T, U parameters (59H62j as well as the 
singlet-induced NLO corrections to the W boson mass as presented in Ref. [ 43 ], 


(3.) perturbativity of the couplings as well as the requirement on the potential to be bounded 
from below, Eqs. Q and ([5]), 

(4.) limits from perturbativity of the couplings as well as vacuum stability up to a certain scale 
H run , where we chose ~ 4 x 10 10 GeV as benchmark point (these constraints will only 
be applied in the high mass region, see Section IIIC for further discussion), 


2 In fact, all Higgs self-couplings depend on tan /3. However, in the factorized leading-order description of production 
and decay followed here, and as long as no experimental data exists which constrains the Higgs boson self-couplings, 
only the Hhh coupling needs to be considered. 
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Ill Theoretical and experimental constraints 


(5.) upper cross section limits at 95% C.L. from null results in Higgs searches at the LEP, Tevatron 
and LHC experiments, 

(6.) consistency with the Higgs boson signal rates measured at the LHC experiments. 

The constraints (1.) - (4.) have already been discussed extensively in a previous publication [ 41 ], 
where the scan was however restricted to the case that mu > 600 GeV. In the following, we will 
therefore briefly recall the definition of the theoretically motivated bounds and comment on their 
importance in the whole mass range m^/u E [1, 1000] GeV. 


A. Perturbative unitarity 


Tree-level perturbative unitarity [631 541 puts a constraint on the Higgs masses via a relation 
on the partial wave amplitudes aj(s ) of all possible 2 —> 2 scattering processes: 


|Re(aj(s))| < 


( 21 ) 


where the partial wave amplitude ao poses the strongest constraint. Following Ref. HU, we consider 
all 2 ^ 2 processes X, X 2 -> Y 1 Y 2 , with (X lt X 2 ), (V, V>) E (W+ W ~, ZZ, hh, hH, HH), and 


impose the condition of Eq. (21) to the eigenvalues of the diagonalized scattering matrix. Note 
that the unitarity constraint based on the consideration of WlWl WlWl scattering alone, 
leading to mu < 700 GeV (as e.g. in Ref. [ 25 ] ), is much loosened when all scattering channels are 
taken into account [58l [65] . 


In general, perturbative unitarity poses an upper limit on tan/3. In the decoupling case, which 
corresponds to sin a —> 0 (1) for the light (heavy) Higgs being SM-like, it is given by [41] 


1 

tan 2 /3 < -—b O (a) for 

3 m z 


ao(hh —> hh) < 0.5, 


( 22 ) 


where h and m refer to the purely singlet Higgs state and its respective mass. 

While in the high mass scenario this bound is always superseded by bounds from perturbativity 


of the couplings, cf. Section III C, in the low mass scenario this poses the strongest theoretical bound 


on tan /3. We exemplarily show the upper limits on tan /3 derived from perturbative unitarity in 
Fig.[l]for the low mass range m E [20,120] for the four values of sin a = 1.0, 0.9,0.5,0.0. The bounds 
on tan/3 are strongest for small values of sin a. However, values too far from the decoupling case 
sin a ~ 1 are highly constrained by Higgs searches at LEP as well as by the LHC signal strength 
measurements of the heavier Higgs at ~ 125GeV, cf. Sections HIE and|IIIF for more details. 


B. Perturbativity of the couplings 

For perturbativity of the couplings, we require that 

|Ai| <4 tt, ie (1,2,3). (23) 
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C RGE evolution of the couplings 



m [GeV] 

FIG. 1: Maximally allowed values for tan/3 in the low mass range, m £ [20,120] GeV, for various values of 
sin a = 1.0, 0.9,0.5, 0.0, from considering only perturbative unitarity. 

At the electroweak scale, these bounds do not pose additional constraints on the parameter space 
after limits from perturbative unitarity have been taken into account. 


C. Renormalization group equation evolution of the couplings 


While perturbativity as well as vacuum stability and the existence of a local minimum at the 
electroweak scale are necessary ingredients for the validity of a parameter point, it is instructive 
to investigate up to which energy scale these requirements remain valid. In particular, we study 
whether the potential is bounded from below and features a local minimum at energy scales above 
the electroweak scale. In order to achieve this, we promote the requirements of Eqs. ©> © , and 

) are evolved according to the one-loop 

) 


(23) to be valid at an arbitrary scale /r run , where A* (/j 
renormalization group equations (RGEs) (see e.g. Ref. 


d_ 

dt 


Ai ~ 10 7 r 2 { 4 ^3 + 12 ~ l ~ 6 Ai dt ~ 3 Dt ~ 2 ^ ^ 5 i) + 16 2 ^ 


1 


± x =r 

dt 2 16 7T 2 

4a 3 = t^A 3 


[A§ + 9 A|] 


dt 


16 7T 2 


6 Ai + 3 A 2 + 2 A 3 + 3 y\ — - (3 g 2 + g{) 


(24) 

(25) 

(26) 


Here we introduced t = 2 log (/i TU nA') as a dimensionless running parameter. The initial condi¬ 


tions at the electroweak scale require that Aj (/r res = v) are given by Eqn. (13). The top Yukawa 
coupling y t as well as the SM gauge couplings g, g\ evolve according to the one-loop SM RGEs, 
cf. Appendix |B| For the decoupling case as well as to cross-check the implementation of the running 
of the gauge couplings we reproduced the results of Ref. m- 
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C RGE evolution of the couplings 


Effects of perturbativity and vacuum stability, n rur) = 2.7 x 10 10 GeV 



FIG. 2: Limits in the (sin a, tan /3) plane for mn = 600 GeV from requiring perturbativity and vacuum 
stability at a scale /x mn = 2.7 x 10 10 GeV using RGE evolution. Taken from Ref. HU- 


As in Ref. HD. we require all RGE-dependent constraints to be valid at a scale which is slightly 
higher than the breakdown scale of the SM, such that the singlet extension of the SM improves 
the stability of the electroweak vacuum. The SM breakdown scale is defined as the scale where 
the potential becomes unbounded from below in the decoupled, SM-like scenario. With the input 
values of rrih = 125.14 GeV and v = 246.22 GeV, a top mass of 173.0 GeV as well as a top-Yukawa 
coupling yt(mt ) = 0.93587 and strong coupling constant a s (mz) = 0.1184, we obtain as a SM 
breakdown scale 3 


Mrun, SM bkdw ~ 2.5 x IQ 10 GeV. 


We therefore chose as a slightly higher test scale the value /i run> sta b ~ 4.0 x 10 10 GeV. Naturally, 
we only apply this test to points in the parameter space which have passed constraints from 
perturbative unitarity as well as perturbativity of the couplings at the electroweak scale. Changing 
the scale to higher (lower) values leads to more (less) constrained regions in the models parameter 
space fib ]. 

In the high mass scenario we see the behavior studied in Ref. m for Higgs masses > 600 GeV 
continuing to the lower mass ranges. The strongest constraints that impact different parts of the 
(sin a, tan /3) parameter space are displayed in Fig. [ 2 ] for a heavy Higgs mass of run = 600 GeV 
(taken from Ref. |4Tj). Two main features can be observed: 

First, the upper value of tan/3 for fixed Higgs masses is determined by requiring perturbativity 
of A 2 as well as perturbative unitarity, cf. Section eh A| Second, the allowed range of the mixing 


3 As has been discussed in e.g. Ref. [67], the scale where Ai = 0 in the decoupling case strongly depends on the 
initial input parameters. However, as we are only interested in the difference of the running in the case of a 
non-decoupled singlet component with respect to the Standard Model, we do not need to determine this scale to 
the utmost precision. For a more thorough discussion of the behavior of the RGE-resulting constraints in case of 
varying input parameters, see e.g. Ref. HD- 
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C RGE evolution of the couplings 


angle sin a is determined by perturbativity of the couplings as well as the requirement of vacuum 
stability, especially when these are required at renormalization scales // run , which are significantly 
larger than the electroweak scale. Small mixings are excluded by the requirements of vacuum 
stability 4 as well as minimization of the scalar potential. This corresponds to the fact that we 
enter an unstable vacuum for fv TVLn > /isM,bkdw for sin a ~ 0. 

In summary, the constraints from RGE evolution of the couplings pose the strongest bounds 
on the minimally allowed value | sina| and the maximal value of tan /3 in the high mass scenario. 
Note that, for lower m#, the (sin a, tan /3) parameter space is less constrained, as will be discussed 
in Section IIV Al 

In the low mass scenario, i.e. where the heavier Higgs state is considered to be the discovered 
Higgs boson, none of the points in our scan fulfilled vacuum stability above the electroweak scale. 
This is due to the fact that for a relatively low nih, the value of Ai at the electroweak scale is quite 
small, cf. Eq. (13). In the non-decoupled case, | sin o:\ 1, Ai then receives negative contributions 

in the RG evolution towards higher scales, leading to Ai (// nm ) < 0 already at relatively low scales 
/i run , corresponding to the breakdown of the electroweak vacuum. Hence, in the low mass scenario, 
the theory breaks down even earlier than in the SM case. In the analysis presented here, we will 
therefore refrain from taking limits from RGE running into account in the low mass scenario. Then, 
the theoretically maximally allowed value of tan/3 is determined from perturbative unitarity and 
rises to quite large values, where we obtain tan/3 max < 50, depending on the value of the light 
Higgs mass rrih ■ 

Further constraints on tan /3 in the low mass scenario stem from the Higgs signal rate observables 
through the potential decay H 


hh, as will be discussed in Section IIIF 


D. The W boson mass and electroweak oblique parameters 5, T, U 

Recently, the one-loop corrections to the W boson mass, mw, for this model have been calcu¬ 
lated in Ref. [23]. In that analysis, mw is required to agree within 2a with the experimental value 
rnff? = 80.385 ± 0.015 GeV RjUffTT] . leading to an allowed range for the purely singlet-induced 
corrections of A m^l s 6 [—5MeV; 55MeV]. Theoretical uncertainties due to contributions at even 
higher orders have been estimated to be O (IMeV). The one-loop corrections are independent 5 of 
tan/3 and give rise to additional constraints on | sina|, which in the high mass scenario turn out to 
be much more stringent [23] than the constraints obtained from the oblique parameters S, T and 

U [59H62]. 


4 For the requirement of vacuum stability, we found that in some cases the coupling strengths vary very mildly over 
large variations of the RGE running scale. In these regions the inclusion of higher order corrections in the spirit 
of Ref. [67] seems indispensable. Therefore, all lower limits on the mixing angle originating from RGE constraints 
need to be viewed in this perspective. In fact, such higher order contributions to the scalar-extended RGEs have 
recently been presented in Ref. [6S]. However, the authors did not specifically investigate the higher order effects 
on parameter points which exhibit small variations over large energy scales at NLO. 

5 In the electroweak gauge sector, tan /3 only enters at the 2-loop level when the Higgs mass sector is renormalized 
in the on-shell scheme. 
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D The W boson mass and S, T, U 



m [GeV] 


FIG. 3: Maximal allowed values for |sina| in the high mass region, m# £ [130,1000] GeV, from NLO 
calculations of the W boson mass {red, solid ) [12, electroweak precision observables (EWPOs) tested via 
the oblique parameters S , T and U {orange, dashed), as well as from the perturbativity requirement of the 
RG-evolved coupling Ai {blue, dotted), evaluated at tan ,5 = 0.1. For Higgs masses mu < 800 GeV the 
NLO corrections to the W boson mass yield the strongest constraint. 


Figure [3] shows the maximally allowed mixing angle obtained from the mw constraint as a 
function of the heavy Higgs mass mn in the high mass scenario. For comparison, we also included 
the limit stemming from the electroweak oblique parameters S, T, and U (see below), as well as 
from requiring perturbativity of Ai (/r m n), evaluated at tan (3 = 0.1. We see that for m# < 800 GeV 
constraints from yield the strongest constraint. The oblique parameters S, T and U do not 
pose additional limits on the allowed parameter space. 

In the low mass region, as discussed in Ref. [32], the NLO contributions within the Higgs singlet 
extension model even tend to decrease the current ~ 20 MeV discrepancy between the theoretical 
value rriw in the SM m and the experimental measurement [SMI]. However, substantial re¬ 
duction of the discrepancy only occurs if the light Higgs has a sizable doublet component. Hence, 
this possibility is strongly constrained by exclusion limits from LEP and/or LHC Higgs searches 
(depending on the light Higgs mass) as well as by the LHC Higgs signal rate measurements. 

In the low mass region the electroweak oblique parameters pose non-negligible constraints, as 
will be shown in Section IV B However, these constraints are again superseded once the Higgs 


signal strength as well as direct search limits from LEP are taken into account, cf. Section HIE 
and III F| respectively. 

In our analysis we test the constraints from the electroweak oblique parameters S, T and U by 
evaluating 

Xstu = x r C _1 x, 

with x 3 ={S — S,T — T,U — U), where the observed parameters are given by [73j 


(27) 


S = 0.05, T = 0.09, U = 0.01, 


(28) 
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D The W boson mass and S, T, U 


and the unhatted quantities denote the model predictions [43j. 6 The covariance matrix reads m 


(C)ii 


/ 0.0121 0.0129 -0.0071 

0.0129 0.0169 -0.0119 

v -0.0071 -0.0119 0.0121 


(29) 


We then require Xstu — 8.025, corresponding to a maximal 2 a deviation given the three degrees 
of freedom. 


E. Exclusion limits from Higgs searches at LEP and LHC 


Null results from Higgs searches at collider experiments limit the signal strength of the second, 
non SM-like Higgs boson. Recall that its signal strength is given by the SM Higgs signal rate scaled 
by (cos a) 2 in the low mass region and, in the absence of Higgs-to-Higgs decays, (sin a) 2 in the 
high mass region. Thus, the exclusion limits can easily be translated into lower or upper limits on 
the mixing angle | sina|, respectively. 7 

We employ HiggsBounds-4.2.0 |52fl55] to derive the exclusion limits from collider searches. 
The exclusion limits from the LHC experiments 8 are usually given at the 95% C.L.. For most of 
the LEP results we employ the y 2 extension [55] of the HiggsBounds package. 9 The obtained y 2 
value will later be added to the y 2 contribution from the Higgs signal rates, cf. Section IIIF to 
construct a global likelihood. 

The 95% C.L. excluded regions of | sina| derived with HiggsBounds are shown in Fig. [d] as a 
function of the second Higgs mass, assuming a vanishing decay width of the Higgs-to-Higgs decay 
mode H —>• hh. Since all Higgs boson production modes are reduced with respect to their SM 
prediction by a universal factor, limits from LHC Higgs search analyses for a SM Higgs boson can 
be applied straight-forwardly m- In particular, the exclusion limits obtained from combinations 
of SM Higgs boson searches with various final states are highly sensitive. However, so far, such 
combinations have only been presented by ATLAS and CMS for the full 7 TeV dataset [83] E4] 
and for a subset of the 8 TeV data |85]. The strongest exclusions are therefore obtained mostly 
from the single search analyses of the full 7/8 TeV dataset, in particular from the channel H —> 
ZZ — > Ai [86i-S8f in the mass region m e [130,150] GeV and for m > 190 GeV, as well as from 


6 The exact one-loop quantities from Ref. [43] render qualitatively the same constraints as the S, T , U values used 
in Ref. [fTj, which were obtained from rescaled SM expressions [74 1. 

7 Here we neglect the possible influence of interference effects in the production of the light and heavy Higgs boson 
and its successive decay. Recent studies [441IT5H82J have shown that interference and finite width effect can lead to 
sizable deviations in the invariant mass spectra of prominent LHC search channels such as gg —> H —> ZZ* —> 4H 
in the high mass region and thus should be taken into account in accurate experimental studies of the singlet 
extended SM at the LHC. However, the inclusion of these effects is beyond the scope of the work presented here. 

lS HiggsBounds also contains limits from the Tevatron experiments. In the singlet extended SM, however, these limits 
are entirely superseded by LHC results. 

9 The LEP \ 2 information is available for Higgs masses > 4 GeV. For lower masses, we take the conventional 
95% C.L. output from HiggsBounds. 
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Ill Theoretical and experimental constraints 


E Exclusion limits from Higgs searches... 



m [GeV] 



200 400 600 800 1000 

m [GeV] 


(a) Mass region m £ [1,100] GeV. 


(b) Mass region m £ [100,1000] GeV, assuming a 
vanishing Higgs-to-Higgs decay mode, F H ^hh = 0. 


FIG. 4: 95% C.L. excluded values of |sina| from LEP and LHC Higgs searches, evaluated with 

HiggsBounds-4. 2 . 0 in the mass regions m £ [1,100] GeV (a) and m £ [100,1000] GeV ( b ). We assume a 
vanishing decay width for the Higgs-to-Higgs decay mode, TH^hh = 0, hence the displayed results in the 
high mass region correspond to the most stringent upper limit on | sina| that can be obtained from current 
LHC Higgs searches. The other Higgs boson mass is set to 125.14 GeV and is indicated by the dashed, 
magenta line in Fig. |4(b)| 


the H -> WW —> ivtv channel [10, [89j [90J in the mass region m £ [160,170] GeV due to the 
irreducible ZZ background in the H —> ZZ —> analyses. For Higgs masses m £ [65,110] GeV 
the only LHC exclusion limits currently available are from the ATLAS search for scalar diphoton 
resonances m • However, these constraints are weaker than the LEP limits from the channel 
► HZ —>• ( bb)Z [92], as can be seen in Fig. 4(b) In the remaining mass regions with 


e 1 e 


m > 110 GeV the CMS limit j85] from the combination of SM Higgs analyses yields the strongest 
constraint. For very low Higgs masses, m < 10 GeV, the LEP constraints come from Higgs pair 
production processes, e + e _ —> hh —> t + t~t + t~ and T + r~bb [92j, as well as from the decay-mode 
independent analysis of e + e _ —> Zh by OPAL [93]. The latter analysis provides limits for Higgs 
masses as low as 1 keV. 

In the presence of Higgs-to-Higgs decays, BR(iL —>• hh) / 0, additional constraints arise. In case 
of very low masses, m < 3.5 GeV, these stem from the CMS search in the H —> hh —> /x + /r _ /r + ^ _ 
channel [[91] . and for large masses, m £ [260,360] GeV, from the CMS search for H —> hh with 
multileptons and photons in the final state 123- These limits will be discussed separately in 
Section IV Note that the limit from SM Higgs boson searches in the mass range m > 250 GeV, as 


presented in Fig. |4(b)[ will diminish in case of non-vanishing BR (H —> hh) due to a suppression of 
the SM Higgs decay modes. We find in the full scan (see Section |TV| ) that, in general, BR(iL —> hh) 
can be as large as ~ 70% in this model. Neglecting the correlation between sin a and BR(iL —> hh) 
for a moment, such large branching fractions could lead to a reduction of the upper limit on | sina| 
obtained from SM Higgs searches by a factor of 1/y / 1 — BR(iL —> hh) < 1.8. However, once 
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Ill Theoretical and experimental constraints 


E Exclusion limits from Higgs searches... 


experiment channel 


obs. signal rate obs. mass [GeV] 

ATLAS 

h -> WW -s- Mu 

m 

1 08 +0 - 22 
±.UO_0.20 

- 

ATLAS 

ZZ ^ 0 


1 44+ 0 - 40 

124.51 ±0.52 

ATLAS 

h —> 77 [J; 


1 - 17 + 0.27 

J-- 1 * -0.27 

125.98 ±0.50 

ATLAS 

h —y tt [ST] 


1 40+0-43 
1 - 4Z —0.37 

- 

ATLAS 

Vh V(bb) [6] 


n cri +0.40 
U ' Oi -0.37 

- 

CMS 

h -» WW -s- Mu 

ms 

n 79+0.20 
u * 1 z -0.18 

- 

CMS 

h , ZZ -> U 0 


0.93i°;i 

125.63 ±0.45 

CMS 

h^ 77 0 


l-14±g;i| 

124.70 ± 0.34 

CMS 

h —>■ TT [8] 


0-78±g;i; 

- 

CMS 

Vh V(bb) 0 


i-oo±8;i8 

- 


TABLE I: Higgs boson signal rate and mass observables from the LHC experiments, as implemented in 
HiggsSignals-1.3.0 and used in this analysis. For the mass measurements we combined the systematic 
and statistical uncertainty in quadrature. 


all other constraints (in particular from the NLO calculation of mw ) are taken into account, only 
BR (H —»• hh ) values of up to 40% are found, see Section IV A Fig. |ll(b)| Moreover, in the mass 
region mn ~ 270 — 360 GeV where the largest values of BR(iL —> hh) appear, the mw constraint on 
| sinct| is typically stronger than the constraints from SM Higgs searches, even if BR (H —> hh) = 0 
is assumed in the latter. Therefore, given the present Higgs search exclusion limits, the signal rate 
reduction currently does not have a visible impact on the viable parameter space 10 . 


F. Higgs boson signal rates measured at the LHC 

The compatibility of the predicted signal rates for the Higgs state at ~ 125 GeV with the latest 
measurements from ATLAS @H3 EH [97] and CMS [THE] is evaluated with HiggsSignals-1.3.0 
by means of a statistical y 2 measure. The implemented observables are listed in Tab. [Ij In the 
following we denote this y 2 value by y^g, which also includes the y 2 contribution from the Higgs 
mass observables evaluated within HiggsSignals. The latter, however, only yields non-trivial 
constraints on the parameter space if the fit allows a varying Higgs mass in the vicinity of 125 GeV. 
In the low mass scenario, where one of the Higgs bosons is within the kinematical range of the LEP 
experiment, the y 2 value obtained from the HiggsBounds LEP y 2 extension, denoted as XleP’ is 
added to the HiggsSignals y 2 to construct the global likelihood 

Xtot = Xhs + Xlep- (3°) 

The 68% and 95% confidence level (C.L.) parameter regions of the model are approximated by the 
y 2 difference to the minimal y 2 value found at the best-fit point, Ay 2 = y 2 — y 2 nin , taking on values 
of 1 (2.30) and 4 (6.18) in the case of a 1 (2)-dimensional projected parameter space, respectively. 


10 Note, however, that this may change in future with significantly improved exclusion limits from SM Higgs searches. 
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F Higgs boson signal rates... 



-1 -0.5 0 0.5 1 


since 

(a) (sina, tan/3) plane. 


FIG. 5: Branching ratio BR (H —> hh ) in the (sin a, 
niH = 125.14 GeV. It becomes minimal for either sin 



-1 -0.95 -0.9 -0.85 -0.8 -0.75 


sina 

(b) Zoomed region of the (sina, tan/3) plane. The gray 
contours indicate the 1, 2 and 3<r regions preferred by 
the signal rates; the green, dashed line displays the 
95% C.L. limit from LEP, cf. Section [Hill 

in /3) plane for fixed Higgs masses = 50 GeV and 
= 0, cos a = 0 or tan/3 = — cos a/ sina. 


When both Higgs masses are fixed, the fit depends on two free parameters, namely sina and 
tan/3. The latter can only influence the signal rates of the Higgs boson H if the additional decay 
mode H —>• hh is accessible. The branching fraction BR (H —> hh) then leads to a decrease of 
all other decay modes and hence to a reduction of the predictions for the measured signal rates, 
cf. Eq. (16). The sensitivity to tan/3 via the signal rate measurements is thus only given if the 
heavier Higgs state is interpreted as the discovered particle, win ~ 125 GeV (low mass region), and 
the second Higgs state is sufficiently light, wih /$ 62 GeV. If the H —>• hh decay is not kinematically 
accessible, or in the case where the light Higgs is considered as the discovered Higgs state at 
~ 125 GeV, there are no relevant experimental constraints on tan /3. 

In the low mass region, the Higgs signal rate measurements constrain the modulus of the mixing 
angle sin a to be close to 1, such that the heavy Higgs boson has nearly the same coupling strengths 
as the SM Higgs boson. Moreover, in order to obtain sizable predictions for the measured signal 
rates, the branching ratio BR (H —> hh) must not be too large. We illustrate its dependence on 
tan/3 in Fig. [5j where we exemplarily show the branching ratio BR (H —> hh) in the (sina, tan/3) 
plane for fixed Higgs boson masses of nrih = 50 GeV and win = 125.14 GeV. As can be seen 
in Fig. 5(a), the decay H —>• hh is dominant over large regions of the parameter space, with the 
exception of three distinct cases: The branching ratio BR (H —> hh) exactly vanishes in the case 
that 


( i) sina = 0, ( ii) cosa = 0 or (Hi) tan/3 = —cos a/sin a. (31) 

In the first case (i) all couplings of the heavy Higgs boson to SM particles vanish completely, thus 
this case is highly excluded by observations. The second case (ii) corresponds to the complete 
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F Higgs boson signal rates... 



-1 -0.95 -0.9 -0.85 -0.8 -0.75 

since 


(a) (sin a, tan (3) plane. The black solid and dashed line 
indicate constant values of k 2 h = 1 and 2, respectively. 


The other contour lines are the same as in Fig. 5(b) 



(b) Total width scale factor k 2 h for the la, 2 a and 3a 
regions favored by the Higgs boson signal rates. 


FIG. 6: Total width scale factor k 2 h = r tot /r toti g M for fixed Higgs masses = 50 GeV and m# = 
125.14 GeV: (a) in the (since, tan/3) plane and (b) for the regions favored by the Higgs boson signal rate 
measurements only. 


decoupling of the lighter Higgs boson, such that the heavier Higgs is identical to the SM Higgs 
boson. In the third and more interesting case (in) the branching fraction can be expanded in 
powers of (tan/3 + cos a/ sinct): 


BR (H hh) = 


i 4m l 

-L -9— 


m 


I 

h ^ 2 


cos 2 a sin 4 a 


cos a 


sina 


+ 0 


tan (3 + 


cos a 


sm a 


8 ttv 2 m H r S M, tot(m H ) 

where Tsm, tot(mjy) is the total width in the SM for a Higgs boson at mass run- 

In Fig. |5(b)| we show a zoom of the (sina, tan/3) plane, focussing on the low-BR(iL —» hh) 
valley and sina values close to —1. We furthermore indicate the parameter regions which are 
allowed at the 1, 2 and 3<r level by the Higgs signal rate measurements by the gray contour lines. 
The maximally values of BR (H —> hh) ~ 26% allowed by the Higgs signal rate measurements at 
95% C.L. are found for sina very close to —1 and large tan/3 values, i.e. in the vicinity of case 
(ii) discussed above. In the given example with rn^ = 50 GeV, the 95% C.L. exclusion from LEP 


searches, as discussed in Section HIE (cf. Fig. 4(a)), imposes sina < —0.985 and is indicated in 
Fig. |5(b)| by the green, dashed line. 

Finally, we plot the total width scaling factor, defined by n 2 H = r t ot/rtot, SM, in the (sina, 


tan/3) plane in Fig. 6(a) We furthermore plot njj for the parameter regions favored by the Higgs 


signal rates in Fig. 6(b) The largest values of k 2 h allowed by both the signal rates and the LEP 


constraints at 95% C.L. are obtained for sina close to —1. In the example of m.h = 50 GeV 
discussed here, the total width is increased by up to around 34% with respect to the SM. This 
maximal value of total width enhancement is independent of the light Higgs mass (assuming that 
the channel H —> hh is kinematically accessible). 
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Ill Theoretical and experimental constraints 


F Higgs boson signal rates... 



80 100 120 140 160 180 

m [GeV] 


FIG. 7: Ay 2 distribution from the Higgs signal rate observables, obtained from HiggsSignals-1.3.0, as a 
function of the second Higgs boson mass to and the mixing angle sin a. The mass of the other Higgs boson 
is fixed at m = 125.14 GeV, indicated by the dashed, magenta line. The gray contour lines indicate the 
favored parameter space at 68.3%, 95.5% and 99.7% C.L., solely based on the Higgs signal rate observables. 
The green striped (orange patterned) region is excluded by LEP (LHC) Higgs searches at 95% C.L., cf. also 
Fig. a For Higgs masses to below 100 GeV and beyond around 152 GeV the signal rate constraints are 
independent of m. 


We now want to draw the attention to the intermediate mass range, where both mass eigenstates 
can contribute to the signal strength measurements at the LHC. If the masses of the two Higgs 
bosons are well separated, the signal yields measured in the LHC Higgs analyses can be assumed 
to be solely due to the one Higgs boson lying in the vicinity of the signal, m ~ 125 GeV. However, 
in analyses with a poor mass resolution, as is typically the case in search analyses for the decay 
modes H — > W + W~, H -4 t + t~ and VH - 4 - bb, the signal contamination from the second 
Higgs boson needs to be taken into account if its mass is not too far away from 125 GeV. While 
a proper treatment of this case can only be done by the experimental analyses, HiggsSignals 
employs a Higgs boson assignment procedure to approximately account for this situation [56]. 
Based on the experimental mass resolution of the analysis and the difference between the predicted 
mass and the mass position where the measurement has been performed, HiggsSignals decides 
whether the signal rates of multiple Higgs states need to be combined. Hence, superpositions 
of the two Higgs signal rates considered here are possible if the second Higgs mass lies within 
100 GeV < m < 150 GeV. 

In Fig. [ 7 ] we show the HiggsSignals Ay 2 value obtained from the signal rate observables as 
a function of the second Higgs boson mass m and the mixing angle sin a. The mass of the other 
Higgs boson is fixed at m = 125.14 GeV. The scan range for m extends over both the low mass 
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IV Results of the Full Parameter Scan 


and high mass region. Since the Higgs boson at ~ 125 GeV needs sufficiently large signal rates 
to accommodate for the observed SM-like Higgs signal strength, small (large) values of sin a are 
favored in the high (low) mass region, such that the second Higgs boson is rather decoupled. 
We furthermore show the parameter space excluded at 95% C.L. by LEP and LHC searches, as 
previously discussed in Section HIE in Fig. |4j 

In the case of nearly mass degenerate Higgs bosons, rrih ~ mu = 125.14 GeV, the sensitivity on 
the mixing angle sin a significantly decreases, as the signal rates of the two Higgs states are always 
superimposed. There remains a slight dependence of the total signal rate on the Higgs masses, 
though, since the production cross sections and branching ratios are mass dependent. Moreover, 
depending on the actual mass splitting and mixing angle, potential effects may possibly be seen 
in the invariant mass distributions of the high-resolution LHC channels pp —> H —> 77 f?j and 
pp —» H —»• ZZ* —>• 4£, at a future linear collider like the ILC (18, 02] or eventually a muon 
collider [ 421 [98]. However, the sensitivity on sin a completely vanishes in the case of exact mass 
degeneracy, rrih = mu, such that the singlet extended SM becomes indistinguishable from the SM. 

The weak Ay 2 dependence on m outside of the mass degenerate region, i.e. for m > 128 GeV 
and m < 122 GeV, is caused by the superposition of the signal rates of both Higgs bosons in some 
of the H — > W + W~ ,t + t~ and bb channels, as discussed above. These structures depend on the 
details of the implementation within HiggsSignals, in particular on the assumed experimental 
resolution for each analysis. For Higgs masses m below 100 GeV and beyond around 152 GeV the 
sin a limit from the signal rates is independent 11 of m. 

We see that for Higgs masses m in the range between ~ 100 GeV and 150 GeV, the constraints 
from the Higgs signal rates are more restrictive than the exclusion limits from Higgs searches at 
LEP and LHC. For lower Higgs masses, m < 100 GeV, the LEP limits (cf. Fig. |4(a) ) generally 
yield stronger constraints on the parameter space. For higher Higgs masses, m £ [150, 500] GeV, 
the direct LHC limits (cf. Fig. |4(b)j ) are slightly stronger than the constraints from the signal rates, 
however, this picture reverses again for Higgs masses beyond 500 GeV, where direct heavy Higgs 
searches become less sensitive. 


IV. RESULTS OF THE FULL PARAMETER SCAN 


In this section we investigate the interplay of all theoretical and experimental constraints dis¬ 
cussed in the previous section on the real singlet extended SM parameter space, specified by 

m = rrih/H, sin a, tan (3. 

We separate the discussion into the high mass, the low mass and the intermediate (or degenerate) 
mass region of the parameter space. In the high and low mass region, we keep one of the Higgs 


11 This statement is only true if the Higgs state at ~ 125 GeV does not decay to the lighter Higgs. As discussed 
above, at low light Higgs masses rrih < mH/2, the branching ratio BR (H —> hh) can reduce the signal rates of the 
heavy Higgs decaying to SM particles. 
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IV Results of the Full Parameter Scan 



m [GeV] 

FIG. 8: Comparison of all constraints on | sin cr| as a function of the heavy Higgs mass m in the high mass 
region. The Ai perturbativity and perturbative unitarity constraint have been evaluated for tan/3 = 0.1. 


masses fixed at 125.14 GeV and vary the other, while in the intermediate mass region we treat both 
Higgs masses as scan parameters. In the following we first present results for fixed mass m in order 
to facilitate the understanding of the respective parameter space in dependence of sin a, tan/3. 
These discussions will then be extended by a more general scan, where all parameters are allowed 
to vary simultaneously. For each of these scans, we generate around 0(10° — 10 6 ) points. We close 
the discussion of each mass region by commenting on the relevant collider phenomenology. 


A. High mass region 


In this section, we explore the parameter space of the high mass region, m G [130,1000] GeV. 
In general, for masses m > 600 GeV, our results agree with those presented in Ref. m- However, 
we obtain stronger bounds on the maximally allowed value of | sin a) due to the constraints from 
the NLO calculation of mw [S], which has not been available for the previous analysis [4T |. As 
has been discussed in Section HID, Fig.[3j the constraints from mw are much more stringent than 
those obtained from the oblique parameters S, T, and U in the high mass region. 

We compile all previously discussed constraints on the maximal mixing angle in Fig. [ 8 j Fur¬ 
thermore, the (one-dimensional) allowed regions in | sina| and tan/3 are given in Tab. [TTJfor fixed 
values of m. 12 Here, the allowed range of | sin a] is evaluated for fixed tan /3 = 0.15 and we explicitly 
specify the relevant constraint that provides in the upper limit on |sina[. We find the following 


12 Note, that the upper limit on | sina| from the Higgs signal rates is based on a two-dimensional Ay 2 profile (for 
floating rrih) in Fig. [8} whereas in Tab. [II] the one-dimensional Ay 2 profile (for fixed roi) is used. This leads to 
small differences in the obtained limit. 
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IV Results of the Full Parameter Scan 


A High mass region 


m [GeV] 

| sin a| 

source upper limit 

(tan /5)max 

1000 

[0.018,0.17] 

Ai perturbativity 

0.23 

900 

[0.022,0.19] 

Ai perturbativity 

0.26 

800 

[0.027,0.21] 

mw at NLO 

0.29 

700 

[0.031,0.21] 

mw at NLO 

0.33 

600 

[0.038,0.23] 

mw at NLO 

0.39 

500 

[0.046,0.24] 

mw at NLO 

0.47 

400 

[0.055,0.27] 

mw at NLO 

0.59 

300 

[0.067,0.31] 

mw at NLO 

0.78 

200 

[0.090,0.43] 

mw at NLO 

1.17 

180 

[0.10,0.46] 

signal rates 

1.30 

160 

[0.11,0.46] 

signal rates 

1.46 

140 

[0.16,0.31] 

signal rates 

1.67 


TABLE II: Allowed ranges for sin a and tan/3 in the high mass region for fixed Higgs masses m. The 
allowed interval of sin a was determined at tan/3 = 0.15. The 95% C.L. limits on sin a from the Higgs 
signal rates are derived from one-dimensional fits and taken at Ay 2 = 4. The lower limit on sin a always 
stems from vacuum stability, and the upper limit on tan/3 always from perturbativity of A 2 , evaluated at 
sin a = 0.1. The source of the most stringent upper limit on sin a is named in the third column. We fixed 
to/j = 125.14 GeV, and the stability and perturbativity were tested at a scale of ~ 4 x 10 10 GeV. 


generic features: For Higgs masses m > 200 — 300 GeV, the W boson mass NLO calculation 
provides the upper limit on | since), at lower masses the LHC constraints at 95% C.L. from direct 
Higgs searches and the signal rate measurements are most relevant. The purely theory-based limits 
from perturbativity of Ai only become important for m > 800 GeV. Furthermore, in the whole 
mass range, the minimal value of | sin a) and the maximal value of tan /3 are determined by vacuum 
stability and perturbativity of the couplings. 

The corresponding (two-dimensional) allowed regions in the (sina, tan/3) plane for fixed Higgs 
masses are shown in Fig. [9j Their shapes are largely dictated by the perturbativity and vacuum 
stability requirements of the RGE evolved couplings, thus basically resembling the features observed 
before in Fig. [2] for run = 600 GeV. Here, however, the maximally allowed values for the mixing 
angle sin a stem now from the NLO calculation of mw or, at rather low masses mn < 200 GeV, 
from the Higgs signal rates and/or exclusion limits. In all cases, the upper limit on tan/3 stems 
from the perturbativity requirement of RGE evolved couplings. For the degenerate case, rrih ~ 
mu ~ 125 GeV, we a priori find no upper or lower limit on the mixing angle. In the degenerate 
case we do not take limits from RGE running into account, hence the only constraint stems from 
perturbative unitarity which renders an upper limit on tan /3. 

We now extend the discussion and treat run as a free model parameter. The results are presented 
as scatter plots using the following color scheme: 

• light gray points include all scan points which are not further classified, 

• dark gray points fulfill constraints from perturbative unitarity, perturbativity of the cou- 
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IV Results of the Full Parameter Scan 


A High mass region 



sina 


sina 


(a) Higgs masses below 200 GeV. 


(b) Higgs masses above 200 GeV. 


FIG. 9: Allowed regions in the (sina, tan (3) plane in the high mass region for fixed Higgs masses to. For 
m H ^ 200 GeV, the upper limit on the mixing angle stems from mw, while for to# < 200 GeV the upper 
limit is given by the signal strength measurements as well as experimental searches (cf. also Fig. [8]). 


plings, RGE running and the W boson mass, as discussed in Sect. |III| A-D, 

• blue points additionally pass the 95% C.L. exclusion limits from Higgs searches, 

• red/ yellow points fulfill all criteria above and furthermore lie within a 1 / 2 a regime favored 
by the Higgs signal rate observables. 


The results are presented in Fig. 10 in terms of two-dimensional scatter plots in the three 
scan parameters. The point distribution in the (sinct, tan/3) plane shown in Fig. |10(a) neatly 
resembles the features of Fig. [9] discussed above: Small mixings are forbidden from the requirement 
of vacuum stability, while the maximal value for the mixing angle, | sina| < 0.50, is limited by the 
Higgs signal rate observables. Fig. 10(b) illustrates how the upper limit on tan /3, which stems from 
the perturbativity requirement of A 2 , roughly follows the expected ~ scaling, cf. Eq. (22). 

Finally, we can easily recognize the upper limit on the mixing angle sina from the mw constraint 
and the perturbativity requirement of Ai, cf. Fig. [3j in the point distribution in the (m#, sina) 
plane shown in Fig. |10(c) These constraints provide the most stringent upper limit on |sina| 
for Higgs masses m# > 260 GeV. At lower Higgs masses, the upper limit is set by the Higgs 
signal rate measurements and exclusion limits from Higgs searches at the LHC, cf. Fig. | 8 ] Here 
it is interesting to see that the favored | sina| region at Higgs masses rrifj between 130 GeV and 
~ 152 GeV is more restricted than at higher Higgs masses. Two effects play a role here: Firstly, 
the lower limit on | sina| from the vacuum stability requirement is stronger than at larger Higgs 
masses; And secondly, the heavy Higgs boson lies still in the vicinity of the discovered Higgs state, 
such that their signal rates are combined in the H —> rr, H —> WW and VH —> V(bb ) channels, 
where the mass resolution is poor. In total, these channels however favor a slightly lower signal 
strength than obtained for a SM Higgs, thus the fit slightly prefers larger Higgs masses m#, where 
the signal rates are not added for these observables within HiggsSignals, cf. also Fig. [7} 

We now turn to the discussion of the collider phenomenology of the high mass region. Exper- 
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IV Results of the Full Parameter Scan 


A High mass region 



sin a 


(a) (sin a, tan/3) plane. 



200 300 400 500 600 700 800 900 1000 
mn [GeV] 

(b) ( rriH , tan/3) plane. 


1.0 



mn [GeV] 

(c) ( rriH , sin a) plane. 


FIG. 10: Two-dimensional parameter correlations between mu , tan/3 and sin a in the high mass region. 
See text for a description of the color coding. 


imentally, the model can be probed by searches for a SM-like Higgs boson with a reduced signal 
rate and total decay widths, or by direct searches for the Higgs-to-Higgs decay mode H —> hh, 
where h is the light Higgs boson at around 125 GeV. 

We show the allowed values of the branching ratio BR (H —> hh), given by Eq. (16), in Fig. 11 


In Fig. 11(a) we show the dependence on sin a exemplarily for fixed Higgs masses run, whereas 
the full mu dependence is displayed in Fig. 11(b) using the same color code as above. We observe 
that the maximal values of BR (H —» hh) are ~ 40%, reached for large, positive sin a values m, 
and low Higgs masses m ~ 300 GeV. At higher Higgs masses the branching ratio BR {H —> hh) is 
around 20% or slightly higher. 

The LHC production cross section of the heavier Higgs boson H is given by the SM Higgs pro- 
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IV Results of the Full Parameter Scan 


A High mass region 




200 


m# [GeV] 


(a) Allowed values of BR (H —> hh) as function of sin a 
for fixed Higgs masses mu- 


(b) BR (H — > hh) as a function of rriH, using the 
same color code as in Fig. [To] 


FIG. 11: Allowed branching ratios of the Higgs-to-Higgs decay channel H —>• hh in the high mass scenario. 



m H [GeV] 


(a) Heavy Higgs signal rate with SM particles in the 
final state. The orange solid (dashed) curves indicate 
the observed (expected) 95% C.L. limits from the 
latest CMS combination of SM Higgs searches |85| . 



mn [GeV] 

(b) Heavy Higgs signal rate with light Higgs bosons in 
the final state. We display the current expected and 
observed 95% C.L. limits from CMS H —> hh searches 
with 77 bb [55] and bbbb [TUP] final states. 


FIG. 12: Collider signal rates of the heavy Higgs boson H decaying into SM particles (a) or into two light 
Higgs bosons, H —> hh , (6), in dependence of the heavy Higgs mass, mn- The rates are normalized to the 
inclusive SM Higgs production cross section at the corresponding mass value [35] 101 ,'. 1 102 ]. 
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IV Results of the Full Parameter Scan 


A High mass region 



rriff [GeV] 


-8 10 1 



m H [GeV] 


(a) Heavy Higgs signal rate with SM particles in the 
final state for the LHC at 8 TeV. The orange solid 
(dashed) curves indicate the observed (expected) 
95% C.L. limits from the latest CMS combination of 
SM Higgs searches [85] . 


(b) Heavy Higgs signal rate with light Higgs bosons in 
the final state for the LHC at 8 TeV. We display the 
current expected and observed 95% C.L. limits from 
CMS H —> hh searches with 7766 m and bbbb ESS] 
final states. 



mn [GeV] 

(c) Heavy Higgs signal rate with SM particles in the 
final state for the LHC at 14 TeV. 



mn [GeV] 

(d) Heavy Higgs signal rate with light Higgs bosons in 
the final state for the LHC at 14 TeV. 


FIG. 13: LHC signal rates of the heavy Higgs boson H decaying into SM particles ( a,c ) or into two light 
Higgs bosons, H —» hh , (b,d), in dependence of the heavy Higgs mass, mn, for CM energies of 8 TeV (a,b) 
and 14 TeV (c,d). 


duction cross section multiplied by sin 2 a. For convenience, we introduce the rate scale factors m 


K = 


a 

CSM 

a 


x BR (H -> SM) = sin 4 a 


K 


• BR (H —> hh) = sin 2 a 


r S M,tot 
Ttot 
2 r H^hh 


(32) 

(33) 


osm r to t 

for the heavy Higgs collider processes leading to SM particles or two light Higgs bosons in the final 
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IV Results of the Full Parameter Scan 


A High mass region 
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(a) Ratio of total width, Rot, and the Higgs mass, 
m H . 


(b) Suppression of the total width with respect the 
total width of a SM Higgs boson at mass mH- 


FIG. 14: Total width, r to t, as a function of the Higgs mass mu- We display the ratio, T tot /mH (a), as well 
as the suppression factor with respect to the SM Higgs width, r to t/rtot, sm (b). We obtain V to t/mH < 0.02, 
as well as a suppression of 25 % or lower of the total width compared to the SM prediction, in agreement 
with Ref. |¥I1| . 


state, respectively. Here, BR (H —> SM) comprises all possible Higgs decay modes to SM particles. 
Note, that k + k' = sin 2 a corresponds to the inclusive heavy Higgs production rate, normalized 
to the inclusive SM Higgs production rate |35[ 1011 102] , 

The predicted signal rates normalized to the SM production cross section, Eqs. (32) and (33), 
are shown as a function of the Higgs mass mu for the high mass region in Fig. 


12 We further¬ 


more display the current 95% C.L. exclusion limits from the latest CMS combination of SM Higgs 
searches [85j, as well as from direct searches for the H —>• hh process with 7766 [99] and bbbb |100| 
final states. We see that at the current stage, the experimental searches with SM-like final states 
yield important constraints for mu £ 300 GeV. As discussed above, at larger masses the upper 
limit on the mixing angle, and thus on the maximal production cross section, stems either from 
mw or from perturbativity. Note, that the displayed CMS limit from the SM Higgs search combi¬ 
nation [55] is only based on < 5.1 fb _1 of 7 TeV and < 12.2 fb -1 of 8 TeV data, thus not exploiting 


the full available data from LHC run 1. Obviously, future LHC searches for a SM-like Higgs boson 
with reduced couplings in the full accessible mass range will play an important role in probing the 
singlet extended SM. The direct searches for the H —> hh process carried out by CMS in the final 
states 77 bb [99] and, in particular, bbbb [100 ] draw near to the allowed region at masses ~ 450 GeV. 
While they do not yield any relevant constraints at the current stage, these searches will become 
important in this model in the upcoming LHC runs, as they are complementary to the SM-like 
Higgs searches. For reference, we also provide the predicted LHC signal cross section for both the 
SM Higgs signatures and H —>• hh signature for CM energies of 8 and 14 TeV in Fig. 13 Note 


that, as discussed earlier (see footnote [7]) , we do not include effects from the interference with the 
Higgs boson at ~ 125 GeV in these predictions. 
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IV Results of the Full Parameter Scan 


A High mass region 


allowed scale factor and total width, u = 4 x 10 10 GeV 

'run 


* 



(a) (T to t, k) plane for fixed Higgs mass values mu- 


(b) (r tc >t, V) plane for fixed Higgs mass values rriH- 



Ttot [GeV] 

(d) (Ttot, V) plane from the full scan. 


Rot [GeV] 

(c) (Tot, A) plane from the full scan. 


FIG. 15: Allowed regions in the (r to t,K) plane ( a,c ) and (r t ot,V) plane ( b,d ). In (a,b) the results are 
shown for various fixed values of mu- a and n' denote the scaling factors for SM-like decays and the new 
physics channel H —> hh, respectively, cf. Eq. (32) and (33). 


In general, the total width of the heavy Higgs boson is of high interest for collider searches. In 
the SM, the width of the SM Higgs boson rapidly rises with its mass. In Ref. m it was shown that 
in the singlet extended SM the total width of the heavy resonance, r to t (nr#), is highly suppressed 
due to the small mixing angle required. The same behavior is observed here. We show the ratio 
Rot/nr#, as well as the suppression of the width, r to t/rsM, in Fig. 14 We see that the total width 
of the heavy Higgs only amounts to up to ~ 20 — 25% at lower masses mu 200 GeV, while it is 
even further suppressed to below 5 — 15% of the SM Higgs width for masses mu > 300 GeV. At 
mu = 1000 GeV, the total width is still below 25 GeV. In comparison to SM Higgs boson of the 
same mass, the total width of these resonances is therefore highly suppressed, which promises to 
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IV Results of the Full Parameter Scan 


A High mass region 


enhance the validity of a narrow width approximation in this mass range . 13 

For completeness, we show the allowed parameter space in the (r to t,/t) and (Ft 0 t, k 1 ) planes 
in Fig. [l5j If these predictions are taken as independent input parameters in future Higgs boson 
collider searches, a direct comparison with the experimental results renders additional constraints 
and — in case of a discovery — could possibly lead to an exclusion of the entire model. 


B. Low mass region 


We now consider the low mass region, i.e. we set the heavy Higgs mass to m# = 125.14 GeV 
and investigate the parameter space with rrih € [1, 120] GeV. In contrast to the high mass region, 
results from LEP searches play an important role in this part of parameter space. As discussed in 
Section [III C] we here do not apply limits from RGE running of the couplings. Before constraints 
from the signal rates are taken into account, this a priori leads to much larger allowed values 
for tan /3, where the upper limit on tan (3 stems from perturbative unitarity, cf. Fig. [l] However, 
whenever the additional decay H —> hh is kinematically allowed, tan (3 values > 1 generally result 


in large branching ratios for this channel, cf. Fig. |5(a) This immediately imposes a quite strong 
suppression of the SM decays of the heavy Higgs state, leading to strong bounds on the minimal 


|sina| value from the signal rates, cf. Section IIIF However, we should keep in mind that in 
parameter regions where tan/3 ~ — cos ct/sin a, the branching ratio for H —>• hh decreases signif¬ 
icantly, thus restoring the signal strength of the heavy Higgs boson to sin 2 a times the SM Higgs 
signal strength. In the mass range where the additional decay is not allowed and up to values of 
m h 100 GeV, the strongest limits on the mixing angle stem from LEP Higgs searches in the 


channel e + e —> Zh —> Z(bb) [92J. For larger Higgs masses, the Higgs signal rates yield stricter 
limits on | sina| than the exclusion limits from LEP and LHC, cf. also Fig.[7J We have summarized 


our finding in the low mass scenario in Tab. Ill 


We now turn to the discussion of the full scan. In order to highlight the importance of LEP 
constraints in the low mass region, we employ the following color coding for the plots: 


• Light gray: points which fail theoretical constraints. 

• Dark gray: points which are excluded by LHC Higgs searches. 

• Blue: points allowed by LHC Higgs searches, but excluded by > 95% C.L. by LEP searches. 

• Dark green: points consistent with LEP constraints within 2a. 

• Light green: points consistent with LEP constraints within lcr. 

• Yellow: points favored within 2er in the global fit (HiggsSignals X 2 + LEP x 2 )- 

• red: points favored within la in the global fit (HiggsSignals X 2 + LEP x 2 )- 


13 See e.g. Ref. [75] for the discussion of finite width effects for SM-likc Higgs bosons in the mass range ran > 200 GeV. 
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IV Results of the Full Parameter Scan 


B Low mass region 


m h [GeV] 

| sill Qi|min, HB 

sina| millj hs 

(tan /3) max 

(tan/3) no H—^hh 

120 

0.410 

0.918 

8.4 

- 

110 

0.819 

0.932 

9.3 

- 

100 

0.852 

0.891 

10.1 

- 

90 

0.901 

- 

11.2 

- 

80 

0.974 

- 

12.6 

- 

70 

0.985 

- 

14.4 

- 

60 

0.978 

0.996 

16.8 

0.21 

50 

0.981 

0.998 

20.2 

0.20 

40 

0.984 

0.998 

25.2 

0.18 

30 

0.988 

0.998 

33.6 

0.16 

20 

0.993 

0.998 

50.4 

0.12 

10 

0.997 

0.998 

100.8 

0.08 


TABLE III: Limits on sin a and tan/3 in the low mass scenario for various light Higgs masses rrih- The 
limits on sin a have been determined at tan/3 = 1. The lower limit on since stemming from exclusion limits 
from LEP or LHC Higgs searches evaluated with HiggsBounds is given in the second column. If the lower 
limit on since obtained from the test against the Higgs signal rates using HiggsSignals results in stricter 
limits, we display them in the third column. The upper limit on tan /3 in the fourth column stems from 
perturbative unitarity for the complete decoupling case (| since| = 1), cf. Fig. [I] In the fifth column we give 
the tan/3 value for which Tn^hh = 0 is obtained, given the maximal mixing angle allowed by the Higgs 
exclusion limits (second column). At this tan/3 value, the |sincn| limit obtained from the Higgs signal rates 
(third column) is abrogated. 


The results are shown in Fig. 16 in terms of two-dimensional scatter plots in the scan parameters. 


In Fig. 16(c), we see that most parameter points allowed by the global fit at the 2a level are found 
for rrih > 80GeV and |sina| > 0.85. For lower Higgs masses the mixing angle is constrained to 
values very close to the decoupling scenario (| sin a) ~ 1). The LEP limits are particularly strong 
in the mass region between 65 GeV and 72 GeV, cf. Fig. |4(a) such that only a few valid points 
are found here, as can be seen best in Fig. |16(b)[ The semi-oval exclusion region in Figs. 16(a) 
and 16(c)| for large tan/3 values and low | sina| values, respectively, corresponds to a 2cr deviation 
in the electroweak oblique parameters S, T and U. 

In Fig. |16(b)[ we observe a drastic change in the distribution of allowed parameter points 
when going to Higgs masses rrih < UT.jy/2 ~ 62 GeV, where the decay mode H —> hh becomes 
kinematically accessible. As discussed earlier in Section III F cf. Fig. [5j the decay H —> hh easily 
becomes the dominant decay mode if tan (3 > 1, unless the mixing angle is very close to | sin a\ = 1. 
Hence, for rrih < mjy/2, most of the allowed points are found for small values of tan/3, since the 
Higgs signal rates favor small values of BR (H —> hh). At larger Higgs masses, rrih > ui-jy/2, the 
favored points are equally distributed over the entire tan /3 range allowed by perturbative unitarity. 


It is interesting to investigate the allowed range of the H —>• hh signal rate in dependence of the 


light Higgs mass. This is shown in Fig. 17(a), where the signal rate is normalized to the SM Higgs 
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IV Results of the Full Parameter Scan 


B Low mass region 



(a) (sina, tan/3) plane. 

1.0 


-0.5 


- 1 . 0 , 



) 40 60 80 

m/, [GeV] 

(b) (mh, tan/3) plane. 



40 60 80 

rn/j [GeV] 

(c) {rrih, sin a) plane. 


FIG. 16: Two-dimensional parameter correlations between to#, tan/? and sin a in the low mass region. 
See text for a description of the color coding. 


boson production. Note, that due to the LEP constraints, the favored points feature a mixing 
angle |sina| ~ 1 and thus the displayed signal rate closely resembles BR (H -» hh ). We see that 
the maximally allowed H —>• hh signal rate is about 22% and is roughly independent on the light 
Higgs mass 14 . This upper limit solely stems from the observed signal rates of the SM-like Higgs 
boson at ~ 125 GeV. These constraints therefore also limit the total width of the heavy Higgs at 
125.14 GeV to values ~ 3 — 5MeV, being in the vicinity of the SM total width of ~ 4.1 MeV. 

We now discuss the case of very low Higgs masses, nth < 4 GeV. Here, the LEP constraints 
stem from the decay-mode independent analysis of e + e - —> Zh by OPAL [93], yielding a slightly 


The reason why the density of allowed points still depends strongly on rrih is that regions which are strongly 
constrained by LEP searches require a large fine-tuning of sin a to render allowed points. 
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B Low mass region 
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(a) Signal rate for the H —> hh signature, normalized 
to the SM Higgs production cross section. 


(b) Signal rate for the H —> hh —> g + g~ g + g~ signature 
at very low masses, normalized to the SM Higgs 
production cross section. The magenta line indicates 
the observed limit from a CMS 8 TeV analysis [94] . 


FIG. 17: Signal rates for the H —► hh signature in dependence of the light Higgs mass: (a) inclusive 
signal rate for Higgs production and successive decay H —► hh , normalized to the SM Higgs production rate; 
(b) exclusive signal rate for the g + g~ g + g~ final state for very low Higgs masses mh- 


weaker limit on the mixing angle, | sina| > 0.965, than at larger masses, cf. Fig. 4(a) In the mass 
region nih £ [1,3] GeV, the branching fraction for the light Higgs decay h —> gg amounts between 
3 — 6%, thus allowing to search for the signature ( pp ) -4 H —» hh —» g + g~g + g~ at the LHC. We 
show the predicted signal rate for this signature 15 for the LHC at a center-of-mass energy of 8 TeV 

, yielding the observed 


in Fig. 17(b) A search for this signature has been performed by CMS 
upper limit 16 displayed as magenta line in the figure. As can be seen, the CMS limit provides 
competitive constraints in this parameter region, excluding a sizable amount of the parameter 
region favored by the global fit. Future LHC searches for the 4 g signature therefore have a good 
discovery potential in this mass region. Other final states, composed of r leptons, strange or charm 
quarks, could be exploited at a future linear e + e~ collider like the ILC. 

A very light Higgs boson h with mass values up to the bb threshold can also be probed at 
R-factories in the radiative decay T —> h'y [103], with successive decay of the light Higgs boson to 
r-lepton, muon or hadron pairs. Here, we provide a rough estimate of the present constraints. 


15 We consider here only the Higgs production via gluon gluon fusion. 

10 This exclusion limit is not provided with HiggsBounds-4.2.0, because the expected limit from the CMS analysis 
is not publicly available. 
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IV Results of the Full Parameter Scan 


B Low mass region 


m h [GeV] 

90% C.L. upper limit on BR(Y(ls) —> / 17 , h —>■...), 
h —> fi + 111 h —v t + t ~ 109 h —> qq [L12 h, —¥ ss 112 

9l 

(upper limit) 

sin a 

(lower limit) 

1.0 

~ 4 

10" 6 


- 


~ 5 

10" 6 



~ 0.25 

~ 0.87 

2.0 

~ 5 

10- 6 


- 


~ 1 

10- 4 

- 5 

10~ 5 

~ 1.16 

- 

3.0 

~ 6 

10- 6 


- 


~ 2 

10- 4 

- 8 

10~ 5 

~ 7.82 

- 

4.0 

~ 8 

10" 6 

1.2 

10 

-5 

~ 4 

10" 4 

- 3 

IO - 4 

2.06 

- 

5.0 

~ 8 

10" 6 

9.1 

10 

-6 

~ 3 

10- 4 

- 3 

10- 4 

0.68 

0.57 

6.0 

~ 1 

10" 5 

2.3 

10 

-5 

~ 5 

10- 5 

- 8 

10~ 5 

1.59 

- 

7.0 

~ 1 

10" 5 

1.6 

10 

-5 

~ 3 

10- 4 

- 1 

10- 4 

1.33 

- 

8.0 

~ 2 

10" 5 

3.2 

10 

-5 

~ 1 

10" 2 

~ 4 

10~ 5 

4.45 

- 


TABLE IV: Constraints from B-factories on a light Higgs boson with mass The second to fifth column 
list the current experimental 90% C.L. upper bounds on the decay rate of T(ls) —> /17 and successive Higgs 
decay (specified in the second title row). The inferred upper limit on the rescaling factor of the bottom 
Yukawa coupling in given in the sixth column, and — if possible — the lower limit on the singlet-doublet 
mixing angle |sina| is given in the last column. We indicate the most relevant constraint for the model 
(yielding the listed limits on the model parameters) by bold numbers. 


The decay rate for the 1 bound state Y(ls) to the Higgs-photon final state (normalized to 
the decay rate of T(ls) —> n + ^~) is given by [103] 


BR(T(ls) ->• hj) _ glG F m 2 b 


BR(T(ls) 




v/2: 


1 - 


ml 


no. 


-^T(la), 


X 7, 


(34) 


where Gp is the Fermi constant, mb the bottom quark mass, a the fine-structure constant and 
BR(Y(ls) —> ~ 2.48 ± 0.05% [ 104] . The factor T represents higher-order corrections. The 

one-loop QCD corrections have been calculated in Ref. [1051 1106] and are known to reduce the 
leading-order estimate by up to 84%, see Ref. HOT) for an extended discussion. In our model, the 
rescaling factor of the bottom Yukawa coupling of the light Higgs is simply given by < 5 ^ = cos a. 

Recent experimental searches have been carried out bv BaBar [108K112] and CLEO fmj, fo¬ 
cussing on the search for a light CP -odd Higgs boson motivated by certain next-to-minimal su¬ 
persymmetric standard model (NMSSM) scenarios [114H116] . The 90% C.L. upper limits on the 
branching fraction of these search signatures are typically of r-*j 0(10 4 — 10 6 ) and are listed for 


representative values of the light Higgs mass in Tab. IV (cf. also Refs. [117H119] for more details). 
Generally, these limits underlie large statistical fluctuations, thus we prefer to use a roughly esti¬ 
mated mean value and indicate this by a in front of the quoted number. Using the SM Higgs 
boson branching ratios for h —> h —> t + t~, h —>• gg and h —> ss in this mass region, we 

can infer a 90% C.L upper limit on the rescaling factor of the bottom Yukawa coupling, < 75 , which 


is listed in Tab. IV If this upper limit is below 1, we furthermore quote the resulting lower limit 
on the mixing angle | sin a| in the table. The resulting limits cannot compete with those obtained 
from direct LEP searches, however, future B-physics facilities such as the Belle II experiment at 
the Super KERB accelerator 1120 will be able to probe the yet unexcluded region. 
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IV Results of the Full Parameter Scan 


B Low mass region 


m h [GeV] 

^ TeV [ P b] 

^ TeV [P b ] 

120 

3.28 

8.40 

110 

3.24 

8.12 

100 

6.12 

14.96 

90 

6.82 

16.26 

80 

2.33 

5.41 

70 

2.97 

6.73 

60 

0.63 

1.38 

50 

0.45 

0.96 

40 

0.74 

1.50 


TABLE V: Maximally allowed cross sections, a gg = (cos 2 «) max x cr SS) SM, for direct light Higgs production 
at the LHC at CM energies of 8 and 14 TeV after all current constraints have been taken into account. The 
SM Higgs production cross sections have been taken from Ref. [351 H21j . 


Finally, we want to comment that despite of the quite strong constraints in the low mass region, 
a substantial number of low mass Higgs bosons could already have been directly produced at the 
LHC. Table [Y] exemplarily lists the maximally allowed LHC cross sections for direct production in 
gluon gluon fusion for a selected range of light Higgs masses at CM energies of 8 and 14 TeV. 1 ' 
We encourage the LHC experiments to explore the feasibility of experimental searches within the 
low mass region and to potentially extend the searches for directly produced scalars into this mass 
range. 


C. Intermediate mass region 

For the intermediate mass region, which contains the special case of mass-degenerate Higgs 
states, we treat both Higgs masses as free parameters in the fit, mh,m,H E [120,130] GeV. Note, 
that the following discussion is based on a few simplifying assumptions about overlapping Higgs 
signals in the experimental analyses. It should be clear that a precise investigation of the near 
mass-degenerate Higgs scenario can only be performed by analyzing the LHC data directly and 
is thus restricted to be done by the experimental collaborations (see e.g. Ref. [7] for such an 
analysis). Nevertheless, we want to point out this interesting possibility here and encourage the 
LHC experiments for further investigations. 

If the Higgs states have very similar masses, their signals cannot be clearly distinguished in the 
experimental analyses and (to first approximation) the sum of the signal rates has to be considered 
for the comparison with the measured rates. Moreover, the observed peak in the invariant mass 
distribution in the H —» 77 and H ZZ* &£ channels, which is fitted to determine the Higgs 
mass, would actually comprise two (partially) overlapping Higgs resonances, where the height of 


17 We thank M. Grazzini for providing us with the production cross sections for rrih < 80 GeV. 
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C Intermediate mass region 


to [GeV] 

sin a hb 

sina| H s 

(tan /3) max 

130 

< 0.806 

< 

0.370 

7.76 

129 

< 0.881 

< 

0.373 

7.81 

128 

< 0.988 

< 

0.377 

7.88 

127 

- 

< 

0.381 

7.94 

126 

- 

< 

0.552 

8.00 

125 

- 


- 

8.07 

124 

- 

> 

0.793 

8.13 

123 

- 

> 

0.864 

8.20 

122 

- 

> 

0.904 

8.26 

121 

- 

> 

0.913 

8.34 

120 

> 0.410 

> 

0.918 

8.41 


TABLE VI: Limits on sin a and tan/3 in the intermediate mass scenario. We fix one Higgs mass at 
125.14 GeV and vary the mass of the other Higgs state, m. The limit on sin a that stems from LHC 
Higgs searches evaluated with HiggsBounds is given in the second column (if available). The limit on sin a 
obtained from the test against the Higgs signal rates with HiggsSignals is given in the third column. Note, 
that depending on the mass hierarchy, we have either an upper or lower limit on sin a, indicated by the 
“<” and “>”, respectively. The upper limit on tan/3 is given in the fourth column and always stems from 
perturbative unitarity, see also Fig. [l] Note, that we do not impose constraints from perturbativity and 
vacuum stability at a high energy scale via RGE evolution of the couplings here. 


each resonance is governed by the corresponding signal strength. Therefore, for each Higgs analysis 
where a mass measurement has been performed, cf. Tab. |TJ we calculate a signal strength weighted 
mean value of the Higgs masses 18 , 


_ fih ' m h + fi H rn H 
m = - 


Hh. + Hh 

to be tested against the measurement, where the SM normalized signal strengths are given by 

Ea e “ a a(m h / H ) X BR a (m h/H ) 


(35) 


Hh/H 


(36) 


Ea ^SM M X BR SM ("0 
Here, rh denotes the mass value hypothesized by the experiment during to signal rate measurement. 
The index a runs over all signal channels, i.e. Higgs production times decay mode, considered in the 
experimental analysis, and e a denotes the corresponding efficiencies. The predicted cross sections a 
and partial widths are obtained from rescaling the respective SM quantities |35j 101L 102] by cos 2 a 
and sin 2 a for h and H , respectively. As mentioned earlier in Section III F the SM normalized 
signal strengths n h / H contain a slight mass dependence 19 since the SM cross sections and branching 
ratios are not constant over the relevant mass range. 


18 Testing overlapping signals of multiple Higgs bosons against mass measurements by employing a mass average 
calculation is the default procedure in HiggsSignals since version 1.3.0. 

19 This mass dependence is neglected per default in HiggsSignals since additional complications arise if theoretical 
mass uncertainties are present. This is however not the case here, since we use the Higgs masses directly as input 
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IV Results of the Full Parameter Scan 


C Intermediate mass region 


We present limits on sin a and tan j3 for various choices of the Higgs mass m in the intermediate 


mass region in Tab. VI We fixed the other Higgs mass to 125.14 GeV. Depending on whether 
or not the Higgs mass m is larger than 125.14 GeV, we obtain either an upper or lower limit on 
sin a from the LHC Higgs search exclusion limits or signal rate measurements, which are listed 
separately. In the case of nearly degenerate Higgs masses, m = 125 GeV, no limit on sin a can 
be obtained, since the Higgs signals completely overlap. We find that no limits from 95% C.L. 
exclusions from Higgs searches can be obtained for Higgs masses within 121 GeV and 127 GeV. 
Moreover, the limits inferred from the signal rates become weaker the closer m is to 125.14 GeV 
due to the signal overlap. In the full intermediate mass region, the limits inferred from the Higgs 
signal rates supersede the limits obtained from null results in LHC Higgs searches. 


The upper limits on tan/3 listed in Tab. VI correspond to the perturbative unitarity bound 
(cf. Fig. [Tj). Similarly as in the low mass region, we do not impose constraints from perturbativity 
and vacuum stability at a high energy scale here. If these were additionally required, tan /3 would 
be limited to values < 1.86 for m > 125.14 GeV. For lower Higgs masses m no valid points would 
be found. It should be noted, however, that the collider phenomenology does not depend on tan/3 
in the intermediate mass region, since Higgs-to-Higgs decays are kinematically not accessible. 


The results from the full four-dimensional scan are presented in Fig. 18 in terms of two¬ 


dimensional scatter plots, using the same color coding as in the high mass region (see e.g. Fig. 10). 


The correlation between the two Higgs masses, Fig. 18(a), shows that allowed parameter points 
with Higgs bosons in the full intermediate mass region are found, however, at least one of the Higgs 
masses is always required to be roughly between 124 GeV and 126.5 GeV. We can furthermore 
learn from Figs. 18(c) and |18(d)| that allowed points with one Higgs mass being below 124 GeV 
(above 126.5 GeV) feature | sin a) values close to 1 (0), such that the other Higgs state at around 
125 GeV has SM Higgs-like signal strengths. In the near-degenerate case, mh ~ mu ~ 125 GeV, 
all mixing angles sin a are allowed and the model appears indistinguishable from the SM at current 
collider experiments. 

Fig. 18 also shows the correlations of tan/3 with the mixing angle sin a, Fig. |18(b)| and the Higgs 
masses, Figs. 18(e) and |18(f)j As stated earlier, tan /3 does not influence the collider phenomenology 
in the intermediate mass range, thus we find allowed parameter points in the full tan/3 range up 
to the maximal value given by perturbative unitarity. 

A direct search for the second Higgs boson in the intermediate mass region at the LHC seems 
challenging. Even if the mass splitting between the two Higgs states is large enough to be re¬ 
solved by the experimental analyses, we expect the second resonance to be much smaller than 
the established signal. Nevertheless we would like to encourage the LHC experiments to perform 
dedicated resonance searches, in particular in the mass region slightly above the current signal, 
mn ~ (125.5 — 126.5) GeV, since in this case larger values of the mixing angle are still allowed 
while an improvement of the vacuum stability at the high-scale may be obtained. More promising 
prospects to resolve the near mass-degenerate Higgs scenario have future experimental facilities 


parameters. The evaluation of the signal strength according to Eq. (361 can be activated in HiggsSignals by 


setting normalize_rates_to_ref erence_position=. True . in the file usefulbits_HS.f90. 
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FIG. 18: Two-dimensional correlations between mh, mn, sin a and tan/3 in the intermediate mass region. 
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V Conclusions 


like the ILC |18| 12] or a muon collider [l2l198) . where the latter provides excellent opportunities 
to measure the mass and the total width of the discovered Higgs boson via a line-shape scan. 


V. CONCLUSIONS 

In this work, we have investigated the theoretical and experimental limits on the parameter 
space of a real singlet extension of the SM Higgs sector, considering mass values of the second 
Higgs boson ranging from 1 GeV to 1 TeV, i.e. within the accessible mass range of past, current 
and future collider experiments. This study complements a previous work m that was restricted 
to mn £ [600 GeV, 1 TeV] and moreover did not include constraints from direct Higgs collider 
searches. In the present work, either the heavy or the light Higgs state can take the role of the 
discovered SM-like Higgs boson at 125 GeV. We found that up to Higgs masses m < 300 GeV, 
exclusion limits from direct Higgs collider searches at LEP and the LHC, as well as the requirement 
of consistency with the measured SM-like Higgs signal rates pose quite strong constraints. At 
higher Higgs masses, strong limits stem from electroweak precision observables, in particular from 
the W boson mass calculated at NLO, as well as from requiring perturbativity of the couplings 
and vacuum stability. The latter two are tested both at the electroweak scale and at a high scale 
/r ~ 4 x 10 10 GeV using the /3-functions of the theory (see e.g. Ref. m and references therein). 

We performed a exhaustive scan in the three model parameters — specified by the Higgs mixing 
angle, the second Higgs mass and the ratio of the Higgs VEVs — and provided a detailed discussion 
of the viable parameter space and the relative importance of the various constraints. We translated 
these results into predictions for collider observables for the second yet undiscovered Higgs boson, 
which are currently investigated by the LHC experiments. In particular, we focussed on the global 
rescaling factor k for the SM Higgs decay modes, the signal rate for the Higgs-to-Higgs decay 
signature H —> hh as well as the total width T of the new scalar. A typical feature of the model 
is that the total width of the new scalar is quite suppressed with respect to the SM Higgs boson 
at such masses. At very light Higgs boson masses below 10 GeV we found that new results from 
LHC searches for the signature H —>• hh —>• 4/r are complementary to LEP Higgs searches and 
thus probe an unexplored parameter region. Also future R-factories should be able to probe these 
parameter regions through the decay T —?• h'y. 

We furthermore investigated the intermediate mass region, where both Higgs masses are between 
120 GeV and 130 GeV, and discussed some of the experimental challenges in probing this scenario. 
Dedicated LHC searches for an additional resonance in the invariant mass spectra of the H —> 77 
(see Ref. [7] for a CMS analysis) and H ZZ* channel in the vicinity of the discovered 

Higgs boson as well as future precision experiments at the ILC or a muon collider may shed more 
light onto this case. 

The discovery of additional Higgs states is one of the main goals of the upcoming runs of the 
LHC. In this model, two distinct and complementary signatures of the second Higgs state arise. 
Firstly, the H —>• hh decay signature, where the best sensitivity for the LHC is obtained for 
heavy Higgs masses between 250 GeV and roughly 500 GeV. These signatures have been recently 
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explored by ATLAS and CMS pJ9j 100. 122 1 but the analyses are not yet sensitive to constrain 
the parameter space. Secondly, Higgs searches designed for a SM Higgs boson are sensitive probes 
of the parameter space. We strongly encourage the experimental collaborations to continue these 
searches in the full accessible mass range. However, some of the features of the second Higgs state 
discussed in this work, such as the strong reduction of the total width, should be taken into account 
in upcoming analyses. Finally, we hope that the predictions of LHC signal cross sections at a CM 
energy of 14 TeV will be found useful for designing some interesting benchmark points for the 
experimental analyses of this model. 
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Appendix A: Minimization and vacuum stability conditions 


In this appendix we briefly guide the reader through the steps from Eq. ([3]) to Eq. Q , using 
the definition of the scalar fields given in Eq. ([0) . We basically follow the discussion as presented 
in Ref. )T11 . 

With the definition of the VEVs according to Eq. 0, the extrema of V are determined using 
the following set of equations: 


^(v,x) = v ■ ( -m 2 + \\v 2 + ^ x 2 ) = 0 


dV 

dv 

dV 

dx 


(v, x) = x ■ I —/i 2 + X 2 X 2 + — v 2 ) = 0 


The physically interesting solutions have v,x > 0: 

A 2 m 2 - 4f^ 2 


v 2 = A 2 . 

A1A2-f- 

2 Ai^ 2 — 4^ m 2 
x = 


AiA 2 - ¥ 


Alternatively, we use Eq. (Al) to eliminate m 2 and fi 2 , leading to 

2 \ 2 . A3 2 2 \ 2 1 ^3 2 

m = Ai v + —x , n = X 2 X + —u . 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 


38 










B RGEs for SM gauge couplings and the top quark Yukawa coupling 


Since the denominator in Eqs. (A3)—{A4) is always positive (assuming that the potential is well- 
defined), the numerators need to be positive as well in order to guarantee a positive-definite non¬ 
vanishing solution for v and x. 

For the determination of the extrema we evaluate the Hessian matrix: 


( d 2 V 


T~L{v, x) 


dv 2 

d 2 V 


\ dvdx 

From this equation, it is straightforward 
Eqs. Q and ([5]) are satisfied. 


d 2 V \ 
dvdx 
d 2 V 
dx 2 / 

to verify that the solutions are minima if and only if 


2\\v 2 A 3 UX 
A3U.T 2A2X 2 


(A6) 


Appendix B: RGEs for SM gauge couplings and the top quark Yukawa coupling 


This section basically follows the discussion in Ref. m- In the SM, all one-loop RGEs for gauge 
couplings are of the form 

dx o 


dt 


= ax . 


The exact analytic solution for this equation is given by 

x(t = t 0 ) 


x(t) = 


where for t = log (p - ) we have 


1 — a x{t = to) (t — to) ’ 


A 


(Bl) 


t - t 0 = 2 log ( — 


For positive values of a, the coupling reaches the Landau pole when the denominator in Eq. (Bl) 
goes to 0; for negative values, x —> 0 for t —> 00 . 


The Yukawa coupling terms are in turn given by 


dx , o 

— = ax + bx 
dt 


with the solution 


_ ^aC'(to)e a{t *°) 

^/l-be 2 a ^- t o)C'{to)' 

. X 2 

with C (to) = ,,° 2 ? where x(t = to) = defines the initial value. For the top quark Yukawa 

CL~\~0 31 q 

coupling we have 
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However, taking the explicit scale-dependence of the SM gauge couplings into account, the above 
solution needs to be modified such that a(t — to) is replaced by f t a(t') dt!. In this work we chose 
to solve the RGE of the top quark Yukawa coupling numerically. 
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